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Abstract 

We reanalyze the issue of generation of the baryon asymmetry at the elec- 
troweak phase transition in the MSSM and compute the baryon asymmetry 
assuming the presence of non-trivial CP-violating phases in the parame- 
ters associated with the left-right stop mixing term and the Higgsino mass 
fi. Making use of the closed time-path (CTP) formalism of nonequilib- 
rium field theory, we write down a set of quantum Boltzmann equations 
describing the local number density asymmetries of the particles involved 
in supersymmetric electroweak baryogenesis. CP-violating sources manifest 
"memory" effects which are typical of the quantum transport theory and are 
not present in the classical approach. Compared to previous estimates, these 
non-Markovian features enhance the final baryon asymmetry by at least two 
orders of magnitude. This means that a CP- violating phase | sin M | as small 
as 1CU 3 (or even smaller) is enough to generate the observed baryon asym- 
metry. 
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1. Introduction and summary 

The presence of unsuppressed baryon number violating processes at high temper- 
atures within the Standard Model (SM) of weak interactions makes the generation of 
the baryon number at the electroweak scale an appealing scenario [p]]. The baryon 
number violating processes also impose severe constraints on models where the baryon 
asymmetry is created at energy scales much higher than the electroweak scale 0. Un- 
fortunately, the electroweak phase transition is too weak in the SM [Q. This means that 
the baryon asymmetry generated during the transition would subsequently be erased 
by unsuppressed sphaleron transitions in the broken phase. The most promising and 
well-motivated framework for electroweak baryogenesis beyond the SM seems to be 
supersymmetry (SUSY). Electroweak baryogenesis in the framework of the Minimal 
Supersymmetric Standard Model (MSSM) has attracted much attention in the past 
years, with particular emphasis on the strength of the phase transition [|J |5|, [| and the 
mechanism of baryon number generation J7l El B Fl- 



Recent analytical [IT], [L2| and lattice computations |L3[ have revealed that the phase 
transition can be sufficiently strongly first order if the ratio of the vacuum expectation 
values of the two neutral Higgses tan (3 is smaller than ~ 4. Moreover, taking into 
account all the experimental bounds as well as those coming from the requirement of 
avoiding dangerous color breaking minima, the lightest Higgs boson should be lighter 
than about 105 GeV, while the right-handed stop mass might be close to the present 
experimental bound and should be smaller than, or of the order of, the top quark mass 

irj. 

Moreover, the MSSM contains additional sources of CP-violation besides the CKM 
matrix phase. These new phases are essential for the generation of the baryon number 
since large CP-violating sources may be locally induced by the passage of the bub- 
ble wall separating the broken from the unbroken phase during the electroweak phase 
transition. Baryogenesis is fuelled when transport properties allow the CP-violating 
charges to efficiently diffuse in front of the advancing bubble wall where anomalous 
electroweak baryon violating processes are not suppressed. The new phases appear in 
the soft supersymmetry breaking parameters associated to the stop mixing angle and 
to the gaugino and neutralino mass matrices; large values of the stop mixing angle 
are, however, strongly restricted in order to preserve a sufficiently strong first order 
electroweak phase transition. Therefore, an acceptable baryon asymmetry from the 
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stop sector may only be generated through a delicate balance between the values of 
the different soft supersymmetry breaking parameters contributing to the stop mixing 
parameter, and their associated CP- violating phases ||. As a result, the contribution 
to the final baryon asymmetry from the stop sector turns out to be negligible. On 
the other hand, charginos and neutralinos may be responsible for the observed baryon 
asymmetry || HI]. Yet, this is true within the MSSM. If the strength of the electroweak 
phase transition is enhanced by the presence of some new degrees of freedom beyond 
the ones contained in the MSSM, e.g. some extra standard model gauge singlets, light 
stops (predominantly the right-handed ones) and charginos/neutralinos are expected to 
give quantitatively the same contribution to the final baryon asymmetry. 

While in the past few years a lot of effort has been devoted to the study of the 
strength of the phase transition within the MSSM, only very recently some attention 
has been paid to the mechanism by which the baryon asymmetry is produced. The 
baryon asymmetry has been usually computed using a number of steps ]7|, [8|, pL4| : 1) 
look for those charges which are approximately conserved in the symmetric phase, so 
that they can efficiently diffuse in front of the bubble where baryon number violation is 
fast, and non-orthogonal to baryon number, so that the generation of a non-zero baryon 
charge is energetically favoured. Charges with these characteristics in the MSSM are the 
axial stop charge and the Higgsino charge, which may be produced from the interactions 
of squarks and charginos and/or neutralinos with the bubble wall, provided a source 
of CP- violation is present in these sectors; 2) compute the CP- violating currents of 
the plasma locally induced by the passage of the bubble wall; 3) write and solve a 
set of coupled differential diffusion equations for the local particle densities, including 
the CP- violating source terms derived from the computation of the current at step 2) 
and the particle number changing reactions. The solution to these equations gives a 
net baryon number which is produced in the symmetric phase and then transmitted 
into the interior of the bubbles of the broken phase, where it is not wiped out if the 
first transition is strong enough. Notice that the CP-violating sources are inserted into 
the diffusion equations by hand only after the CP-violating currents have been defined 
and computed. This procedure introduces some degree of arbitrariness and -indeed- 
different CP-violating sources have been adopted for the stop and the Higgsino sectors 
in the literature J7|, [J. This is not an academic question. Adopting different sources 
leads to different numerical results for the final baryon asymmetry, especially if the 
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sources are expressed in terms of a different number of derivatives of the Higgs bubble 
wall profile and, therefore, in terms of different powers of the bubble wall velocity v u 
and bubble wall width L u . 

This crucial issue has been recently investigated in where it was shown that 
non-equilibrium Quantum Field Theory provides us with the necessary tools to write 
down a set of quantum Boltzmann equations (QBE's) describing the local particle den- 
sities and automatically incorporating the CP-violating sources. The most appropriate 
extension of the field theory to deal with these issues is to generalize the time contour 
of integration to a closed time-path (CTP). The CTP formalism is a powerful Green's 
function formulation for describing non-equilibrium phenomena in field theory, it leads 
to a complete non-equilibrium quantum kinetic theory approach and to a rigorous com- 
putation of the CP-violating sources for the stop and the Higgsino numbers. In this 
way, the latter have been rigorously defined and the level of arbitrariness of the pre- 
vious treatments has been dismissed ||15|| . What is more relevant, though, is that the 
CP-violating sources- and more generally the particle number changing interactions- 
built up from the CTP formalism are characterized by "memory" effects which are typ- 
ical of the quantum transport theory |T6| , 17] . CP- violating sources are built up when 
right-handed stops and Higgsinos scatter off the advancing Higgs bubble wall and CP 
is violated at the vertices of interactions. In the classical kinetic theory the "scattering 
term" does not include any integral over the past history of the system. This is equiv- 
alent to assuming that any collision in the plasma does not depend upon the previous 
ones. On the contrary, the quantum approach reveals that the CP-violating source is 
manifestly non-Markovian. As we shall see, these memory effects enhance the value of 
of the final baryon asymmetry by at least two orders of magnitude with respect to the 
previous results. This means that the lower bound on the CP-violating phases from 
requiring sueccessful baryogenesis is considerably relaxed -phases as large as 10~ 3 or 
smaller are enough to generate the observed baryon symmetry. This has important 
implications for the supersymmetric CP-problem. We will also investigate the struc- 
ture of the kinetic QBE's derived with the CTP formalism. These equations have an 
obvious interpretation in terms of gain and loss processes. However, the equations are 
manifestly non-Markovian and only the assumption that the relaxation time scale of 
the particle asymmetry is much longer than the time scale of the non-local kernels leads 
to a Markovian description. Further approximations lead to the familiar Boltzmann 
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equations. 

The paper is organized as follows. In section 2 we give a brief description of the 
basic features of the non-equilibrium quantum field theory and the CTP formalism. In 
sections 3 and 4 we compute the quantum transport equations for local particle asym- 
metries in the bosonic and fermion case, respectively. Section 5 contains the explicit 
computation of the CP-violating source for the right-handed stop and the discussion 
about how to go from general QBE's to diffusion/Boltzmann equations. Section 6 is 
devoted to the same issues, but for the Higgsino number. Section 7 is devoted to the 
computation of the final baryon asymmetry and comparison to previous results. We 
conclude with an outlook of our findings in section 8. 

2. The Schwinger-Keldysh formalism for non-equilibrium quantum field 
theory 

In this section we will briefly present some of the basic features of the non-equilibrium 
quantum field theory based on the Schwinger-Keldysh formulation [|18|]. The interested 
reader is referred to the excellent review by Chou et al. ]19| for a more exhaustive 
discussion. 

Since we need the temporal evolution of the particle asymmetries with definite initial 
conditions and not simply the transition amplitude of particle reactions, the ordinary 
equilibrium quantum field theory at finite temperature is not the appropriate tool. The 
most appropriate extension of the field theory to deal with nonequilibrium phenomena 
amounts to generalize the time contour of integration to a closed-time path. More 
precisely, the time integration contour is deformed to run from — oo to +oo and back 
to — oo. 

The CTP formalism (often dubbed as in-in formalism) is a powerful Green's function 
formulation for describing non-equilibrium phenomena in field theory. It allows to 
describe phase-transition phenomena and to obtain a self-consistent set of quantum 
Boltzmann equations. The formalism yields various quantum averages of operators 
evaluated in the in-state without specifying the out-state. On the contrary, the ordinary 
quantum field theory (often dubbed as in-out formalism) yields quantum averages of 
the operators evaluated with an in-state at one end and an out-state at the other. 

Because of the time contour deformation, the partition function in the in-in formal- 
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ism for a complex scalar field is defined to be 



J, J 1 " 



G 

/ I ( J+<p-\ 



J. 



Tr |T (^exp |i / (J0 + J f f 
= Tr T + ^exp 
x T_ (exp -i J + Jl^l) 



P 



where the suffic C in the integral denotes that the time integration contour runs from 
minus infinity to plus infinity and then back to minus infinity again. The symbol 
p represents the initial density matrix and the fields are in the Heisenberg picture 
and defined on this closed time contour. As with the Euclidean time formulation, 
scalar (fermionic) fields <fi are still periodic (anti-periodic) in time, but with <p(t,x) = 
4>(t — i(3,x), (3 = 1/T. The temperature appears due to boundary condition, but time 
is now explicitly present in the integration contour. 

We must now identify field variables with arguments on the positive or negative 
directional branches of the time path. This doubling of field variables leads to six 



different real-time propagators on the contour These six propagators are not 

independent, but using all of them simplifies the notation. For a generic bosonic charged 
scalar field <fi they are defined as 



{x, y) 

G> 0, y) 
G%x,y) 
Gl{x,y) 
Gl(x,y) 



-i{(j>{x)(j>\y)), 
-i{<j>\y)4>{x)) , 

0( x > y) G $(x, y) + 6(y, x)G$(x, y), 
6{y,x)G>(x,y) + 9(x,y)G<(x,y), 
G\ — G^ = G^ — G^, Gf(x, y) ■■ 



Lfj, — Uj: — (jj, — <_T„ 



(2) 



where the last two Green functions are the retarded and advanced Green functions 
respectively and 8(x,y) = 9(t x — t y ) is the step function. For a generic fermion field ijj 
the six different propagators are analogously defined as 



G> (x, y) 
(x,y) 
Gl,(x,y) 
Gl(x,y) 
G r ^{x,y) 



-i(i)(x)4)(y)), 
+i(ip(y)*P(x)}, 

9(x, y)G^(x, y) + 6(y, x)G^{x, y), 
6{y,x)G>{x,y) + 8{x,y)G<{x,y), 
G\ — G^ = G^ — G^p, G^(x, y) - 



(3) 



For equilibrium phenomena, the brackets (• • •) imply a thermodynamic average over all 
the possible states of the system. While for homogeneous systems in equilibrium, the 
Green functions depend only upon the difference of their arguments (x, y) — (x — y) 
and there is no dependence upon (x + y), for systems out of equilibrium, the definitions 
(0) and (|3|) have a different meaning. The concept of thermodynamic averaging is now 
ill-defined. Instead, the bracket means the need to average over all the available states 
of the system for the non-equilibrium distributions. Furthermore, the arguments of the 
Green functions (x,y) are not usually given as the difference (x — y). For example, 
non-equilibrium could be caused by transients which make the Green functions depend 
upon (t x ,t y ) rather than (t x — t y ). 

For interacting systems whether in equilibrium or not, one must define and calculate 
self-energy functions. Again, there are six of them: £*, £*, S < , S > , Y7 and S a . The 
same relationships exist among them as for the Green functions in (|2]) and (0), such as 

S r = £' - £< = £> - £*, S a = E* - £> = £< - £*. (4) 

The self-energies are incorporated into the Green functions through the use of Dyson's 
equations. A useful notation may be introduced which expresses four of the six Green 



functions as the elements of two-by-two matrices [20 



G ~ ( G> -G l ) ' S ~ ( E> -£*" ) ' (5) 
where the upper signs refer to bosonic case and the lower signs to fermionic case. For 
systems either in equilibrium or non-equilibrium, Dyson's equation is most easily ex- 
pressed by using the matrix notation 

G(x,y) = G°(x,y) + J d 4 x 3 J d A x 4: G°(x, x 3 )E(x 3 , x 4 )(5(x4, y), (6) 

where the superscript "0" on the Green functions means to use those for noninteracting 
system. This equation appears quite formidable; however, some simple expressions may 
be obtained for the respective Green functions. It is useful to notice that Dyson's 
equation can be written in an alternate form, instead of (|6|), with G° on the right in the 
interaction terms, 

G(x, y) = G°(x, y)+J d 4 x 3 J d A x, G(x, x 3 )i:(x 3 , x,)G°(x 4 , y). (7) 

Equations. (|6|) and (0) are the starting points to derive the quantum Boltzmann equa- 
tions describing the temporal evolution of the CP- violating particle density asymmetries. 
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3. QBE for bosonic particle density asymmetry 



From now on we will adopt the general method of deriving the QBE's provided by 
Kadanoff and Baym JH]]. In this section our goal is to find the QBE for the generic 
bosonic CP-violating current 

( J +( x )) = i(4>\x) <j>{x)) = [n^x), X/,(x) 

The zero-component of this current represents the number density of particles mi- 
nus the number density of antiparticles and is therefore the quantity which enters the 
diffusion equations of supersymmetric electroweak baryogenesis. 

Since the CP-violating current can be expressed in terms of the Green function 

G^(x,y) as 

<W> = - {^-d^)G<(x,y)\ (9) 

the problem is reduced to find the QBE for the interacting Green function G^(x, y) when 
the system is not in equilibrium. This equation can be found from (^) by operating by 
(\J X +m 2 ) on both sides of the equation. Here m represents the bare mass term of the 
field 0. On the right-hand side, this operator acts only on G'l 



(p x +m 2 ) Gf(x, y) = <5 (4) (x, y)h + / rf 4 x 3 S (x, x 3 )G (f> (x 3 , y) 



(10) 



where I is the identity matrix. It is useful to also have an equation of motion for the 
other variable y. This is obtained from (|7|) by operating by (x\ y +m 2 ) on both sides of 
the equation. We obtain 

G<t>{x,y) (p y + m 2 ) = 5 ( - 4 \x,y)T 4 , + J d^G^x, x 3 )S^(x 3 , y). (11) 

The two equations ( ]TD| ) and ( |TT| ) ar e the starting point for the derivation of the QBE for 
the particle asymmetries. Let us extract from fllPf ) and flTll ) the equations of motions 
for the Green function G^(x,y) 



a x + m 2 ) G<(x,y) 



d A x 3 T, t (j) (x,x 3 )G^(x 3 ,y) - Ej(ar, x 3 )Gl(x 3 , y) 



12) 



G$(x,y) (G„ +m< 



d 4 x 3 G^(x,x 3 )Y^(x 3 ,y)-G^(x,x 3 )Y? 4> (x 3 ,y) . (13) 



If we now substract the two equations and make the identification x = y, the left-hand 
side is given by 

drjX) (dn 4 



(a^-af;) G<(x, y ) 



x=y 



dT 



V -J4> , 



(14) 



and the QBE for the particle density asymmetry is therefore obtained to be 
dn^X) 



dT 



+ V -MX) = - / cf £3 



(15) 



x=y 



where we have defined the center-of-mass coordinate system 



X = (T,X) = -(x + y), (t,r)=x-y. 



(16) 



Notice that T now means the center-of-mass time and not temperature. The identifica- 
tion x = y in Eq. flllf) is therefore equivalent to require t — f — 0. 

In order to examine the "scattering term" on the right-hand side of Eq. fll5f) ) the 
first step is to restore all the variable arguments. Setting x = y in the original notation 
of E^(x, x 3 )G c f,(x 3 , y) gives (X, x 3 )(x 3 ,X) for the pair of arguments 

dn^X) 



or 



+ v -MX) 



d*x 3 m(X,x 3 )G<(x 3 ,X) - X<(X,x 3 )G%(x 3 ,X) 
Gt(X,x 3 )X<(x 3 ,X)-G<(X,x 3 )rt(x 3 ,X)] . (17) 



The next step is to employ the definitions in (0) to express the time-ordered functions 
Gi, G^, S^, and in terms of G^, G^, and G^. Then the time integrals are 
separated into whether t 3 > T or t 3 < T and the right-hand side of Eq. ([17|) reads 



d 4 x 3 [6(T - t 3 



yi>/^f< 1 /^~f<yi> yi</^f> /o>yi< 



+ e{t 3 -T) 



v</^< 1 r i < s r< v</^ y < r^<v< 



The term with t 3 > T all cancel, leaving T > £3. Rearranging these terms gives [[15 
dn^X) 



+ V -J^X) = - / d 6 x 3 



dU 



£>(X,2 3 )G<(X3,X) 



dT 

- G>(X,x 3 )S<(x3,X) + G<(X,x 3 )S>(x3,X)-S<(X,x 3 )G>(x3,X)]. (19) 

This equation is the QBE for the particle density asymmetry and it can be explic- 
itly checked that, in the particular case in which interactions conserve the number of 
particles and the latter are neither created nor destroyed, the number asymmetry 
is conserved and obeys the equation of continuity dn ( j ) /dT+ V -j<j> = 0. During the 
production of the baryon asymmetry, however, particle asymmetries are not conserved. 
This occurs because the interactions themselves do not conserve the particle number 
asymmetries and there is some source of CP-violation in the system. The right-hand 
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side of Eq. (Jl9|), through the general form of the self-energy E^, contains all the infor- 
mation necessary to describe the temporal evolution of the particle density asymmetries: 
particle number changing reactions and CP-violating source terms, which will pop out 
from the corresponding self-energy Sep. If the interactions of the system do not violate 
CP, there will be no CP-violating sources and the final baryon asymmetry produced 
during supersymmetric baryogenesis will be vanishing. 

As we shall see, the kinetic Eq. ( |T§D has an obvious interpretation in terms of gain 
and loss processes. What is unusual, however, is the presence of the integral over the 
time: the equation is manifestly non-Markovian. Only the assumption that the relax- 
ation time scale of the particle asymmetry is much longer than the time scale of the 
non-local kernels leads to a Markovian description. A further approximation, i.e. tak- 
ing the upper limit of the time integral to T — > oo, leads to the familiar Boltzmann 
equation. The physical interpretation of the integral over the past history of the sys- 
tem is straightforward: it leads to the typical "memory" effects which are observed in 
quantum transport theory |16 , [L7]] . In the classical kinetic theory the "scattering term" 
does not include any integral over the past history of the system which is equivalent to 
assume that any collision in the plasma does not depend upon the previous ones. On the 
contrary, quantum distributions posses strong memory effects and the thermalization 
rate obtained from the quantum transport theory may be substantially longer than the 
one obtained from the classical kinetic theory. As shown in ]T5| |, memory effects play a 
fundamental role in the determination of the CP-violating sources which fuel baryoge- 
nesis when transport properties allow the CP- violating charges to diffuse in front of the 
bubble wall separating the broken from the unbroken phase at the electroweak phase 
transition. 

Notice that so far we have not made any approximation and the computation is 
therefore valid for all shapes and sizes of the bubble wall expanding in the thermal bath 
during a first-order electroweak phase transition. 
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4. QBE for fermionic particle density asymmetry 



The generic fermionic CP- violating current reads 



( J ip( x )) = (^(xh^(x)} = [n^(x),4(x) 



(20) 



where indicates a Dirac fermion and 7 M represent the usual Dirac matrices. Again, the 
zero-component of this current represents the number density of particles minus the 
number density of antiparticles and is therefore the relevant quantity for the diffusion 
equations of supersymmetric electroweak baryogenesis. 

We want to find a couple of equations of motion for the interacting fermionic Green 
function G^(x,y) when the system is not in equilibrium. Such equations may be found 
by applying the operators \ i jd x —M^j and (i /) y +M^j on both sides of Eqs. @ and 
(0), respectively. Here M represents the bare mass term of the fermion if). We find 



i P x -Mj G^(x,y) = 5 {4) (x,y)h + J d 4 x 3 T li} (x,x 3 )G i} (x 3 ,y), (21) 
G^x,y)(ifi y +M) = -5^(x,y)h- f d 4 x 3 G4x,x 3 )i:4x 3 ,y). (22) 



We can now take the trace over the spinorial indeces of both sides of the equations, sum 
up the two equations above and finally extract the equation of motion for the Green 



function G^ 



Tr 



i fix +i A 



G^{x,y) 



(fx 3 Tr E^(x,x a )G^(x s ,y) - T^(x,x 3 )G^(x 3 ,y) 



- G^(x,x 3 )12 t Jx 3 ,y) + GUx^x^H^x^y) 



(23) 



Making use of the center-of-mass coordinate system, we can work out the left-hand side 
ofEq. (H) 



Tr 



i jd x G> (T, X, t, r) + G> (T, X, t, r)i 



t=f=o 



--^mx)YHX)) 

d 



t=r=0 



Jt 1 

dX^ ^ ' 



(24) 



The next step is to employ the definitions in ([3D to express the time-ordered functions 
G\, Glp, EL and in terms of G^, £^ and G^. The computation goes along the 
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same lines as the analysis made in the previous section and we get [13 



?H^1+^.J^(X) = J d 3 x 3 J^dt 3 Tr[Z>(X,x 3 )G<(x 3 ,X) 

- g>{x,x 3 )i:<(x 3 ,x) + g<(x,x 3 )e>(x 3 ,x)-j:<{x,x 3 )g>(x 3 ,x)]. (25) 

This is the "diffusion" equation describing the temporal evolution of a generic fermionic 
number asymmetry n^. As for the bosonic case, all the information regarding particle 
number violating interactions and CP- violating sources are stored in the self-energy S^. 

5. The QBE for the right-handed stop number 

As we mentioned in the introduction, a strongly first order electroweak phase transi- 
tion can be achieved in the presence of a top squark lighter than the top quark [JTT], [12 



In order to naturally suppress its contribution to the parameter Ap and hence pre- 
serve a good agreement with the precision measurements at LEP, it should be mainly 
right-handed. This can be achieved if the left-handed stop soft supersymmetry breaking 
mass mjr is much larger than Mz- Under this assumption, only the right-handed stops 
contribute to the axial stop charge. The right-handed stop current J~ associated to 
the right-handed stop i# is given by 

Jf R = i (t* R # t H ) • (26) 

The self-energy of the right-handed stop contains the information about all the 
sources which are responsible for changing nj R in the plasma: scattering processes in- 
volving the top quark Yukawa coupling, axial top number violation processes, strong 
sphaleron interactions and the CP-violating source induced by the presence of CP- 
violating phases in the interactions of the right-handed stop with the Higgs background. 

Eq. (|i~9l) is the QBE describing the right-handed stop number asymmetry. Solving 
this equation represents an Herculean task since it is integral and nonlinear. This 
happens because the self-energy functions S > and are also functions of the full 
non-equilibrium Green functions of other degrees of freedom of the system. We can 
make some progress, though. Since we know that nj identically vanishes if there is no 
CP- violating source in the QBE, i.e. in absence of the Higgs configuration describing 
the bubble wall profile, we first decompose the generic Green function G and self-energy 
function £ as 

G = G° + 5G, S = S° + 5S, (27) 
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where G° and S° represent the fully interacting equilibrium Green functions and self- 
energy functions in the unbroken phase (that is in the absence of the Higgs profile). This 
approximation amounts to retaining only the first-order linear response in the "Higgs 
insertion expansion" around the symmetric phase (H^(x)) = Vi{x) = (i = 1,2). This 
is certainly a good approximation for the case under investigation. We will return to 
this point later. 

If we now linearize Eqs. (|5|) and (|7|) and repeat the procedure described in section 
3, we obtain 



dT 



tR t R tR t R t R tR 

tR t R t R t R t R t R 



;28) 



where we have used the fact that n?~ — j~ =0 and therefore the right-handed stop 

tR tR 



t R J t! 

number asymmetry and current are associated to the Green function 5G~ . 

In the spirit of the diagrammatic approach in the Higgs insertion expansion, we can 
now write the shift in the self-energy as 

5Sr = 5£f p + <5Ef nt + • • • , (29) 

tR t R t R V / 

where <5S~ P is the part of the self-energy responsible for the appearance of the CP- 
violating source and accounts the interactions which change the small particle 

tR 

number asymmetry of the right-handed stop originated at the passage of the Higgs 
profile through a given region of space. Here the dots represent other terms in the 
self-energy describing the interactions which do not change the small particle num- 
ber asymmetry of the right-handed stop (elastic scatterings) - they will not give any 
contribution to the right-hand side of the QBE. 

Eq. (p8~D contains in the right-hand side generic expressions like 

5T£ p G^ + S^ nt G^ + £fi 5G 7 (30) 

t R t R t R t R t R tR V I 

and 

G£ 5Sf p + Gfi 5Sf nt + 5G T . (31) 

tR t R t R t R tR tR \ I 

Let us first analyze the CP-violating source for the right-handed stop number induced 

by a novanishing 5S5 P . 

J ° t R 
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-The CP-violating source for the right-handed stop number- 



The interaction which is responsible for the CP-violating source right-handed stop 
Ir, the left-handed stop t L and the two neutral Higgses Hf 2 is given by 

C = h t t L [A t H° 2 - fi*H°) t* R + h.c. (32) 

Here the soft trilinear term A t and the supersymmeric mass term \i are meant to be 
complex parameters so that lm(A t /j,) is nonvanishing. 

At the lowest level of perturbation, the interactions (p5^) induce a contribution to 
the self-energy of the form [|T5|] 

5^> > (x,y)=g cp (x,y)Gl(x,y), 5£ CP <<(x, y) = g CP (x, y)G$(x, y). (33) 

Here 

9cp(x,y) = h\ [A*v 2 (x) - ixv x {x)\ [A t v 2 (y) - ^*v x {y)\ (34) 

and we can safely approximate the exact left-handed stop Green functions Gj l with 
the corresponding fully interacting equilibrium Green functions in the unbroken phase 
G~ . This is because any departure from thermal equilibrium distribution functions is 
caused at a given point by the passage of the wall and, therefore, is 0(v w ). Since we will 
show that the source is already linear in v w , working with thermal equilibrium Green 
functions in the unbroken phase amounts to ignoring terms of higher order in v u . This 
is as accurate as the bubble wall is moving slowly in the plasma. 

If we now insert the expressions (|33| ) and ( |3"4"D into the QBE (p8|), we get the right- 
handed stop CP-violating source |TH 



= -2i J d 3 x 3 J dt 3 [g C p(X,x 3 ) - g C p(x 3 ,X)} 



x Im 



G£' > (X,x 3 )Gj < (x 3 ,X) 



+ 



Ah'f I d A x 3 [ dt 3 1m(Atfi)[v 2 (X)v 1 (x! i )-V2(x3)v 1 (X)] 



x Im 



G^iX.x^G^ix^X) . (35) 



&t vanishes if the relative phase of A t fi is zero and if the ratio tan/3(x) = V2{x)/vi(x) 
is a constant in the interior of the bubble wall. 

Notice that the source is built up integrating over all the history of the system. 
It is exactly this "memory effect" that is responsible for the enhancement of the final 
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baryon asymmetry. Furthermore, the source is constructed when the right-handed stops 
pass across the wall, they first scatter off the wall and are transformed into left-handed 
stops; the latter subsequently suffer another scattering off the wall and are converted 
again into right-handed stops. If CP-violation is taking place in both interactions, a 
nonvanishing CP-violating source pops out from thermal bath. 

In order to deal with analytic expressions, we can work out the thick wall limit and 
simplify the expressions obtained above by performing a derivative expansion 



OO 



= E k\ d^r v ^ x) ^ ~ xr ■ (36) 

The term with no derivatives vanishes in the expansion (|36|), V2(X)vi(X)—vx(X)v2(X) = 
0, which means that the static term in the derivative expansion does not contribute to 
the source . For a smooth Higgs profile, the derivatives with respect to the time 
coordinate and n > 1 are associated with higher powers of v w /L u , where v u and 
are the velocity and the width of the bubble wall, respectively. Since the typical time 
scale of the processes giving rise to the source is given by the thermalization time of the 
stops 1/IY, the approximation is good for values of L^Tj/v^ ^> 1. In other words, this 
expansion is valid only when the mean free path of the stops in the plasma is smaller 
than the scale of variation of the Higgs background determined by the wall thickness, 
L u , and the wall velocity u w . A detailed computation of the thermalization rate of the 
right-handed stop from the imaginary part of the two-point Green function has been 
recently performed in |22|] by making use of improved propagators and including resum- 
mation of hard thermal loopsQ. The thermalization rate has been computed exactly at 
the one-loop level in the high temperature approximation as a function of the plasma 
right-handed stop mass mq R (T) and an estimate for the magnitude of the two- loop con- 
tributions which dominate the rate for small mqr (T) was also given. If mj R (T) > T, 
the thermalization is dictated by the one-loop thermal decay rate which can be larger 
than T 

The term corresponding to n = 1 in the expansion ( PUD gives a contribution to the 
source proportional to the function 

vx(X)d x v 2 (X) - v 2 (X)d x v 1 (X) = v 2 {X)d x (3{X) 1 (37) 



1 The left-handed stop is usually considered to be much heavier than T and its decay width corre- 
sponds to the one in the present vacuum. 

2 For smaller values of mqr (T), when the thermalization is dominated by two-loop effects (i.e. scat- 
tering), Tj r may be as large as 10~ 3 T 
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which should vanish smoothly for values of X outside the bubble wall. Here we have 
denoted v 2 = v \ + v\. Since the variation of the Higgs fields is due to the expansion 
of the bubble wall through the thermal bath, the source will be linear in v u . The 
corresponding contribution tot he CP-violating source reads 



S ¥R (X)=h 2 t lm(A tf ,)v 2 (X)$(X)l 1 



(3* 



where $(X) = d/3(X)/dt x , 



1: 



tR 



dk 



k 2 



2n 2 ujj ujt 
2Re(/i 



tL t R 



L 



+ (l + 2Re(^J)J(a^,r. 



2(lm(4)+Im(^)) G{^ 



R' *l' t L ' 



(39) 



alld f? B (L) = V H + iT TbW 

tions in the unbroken phase?]. 

The functions J and G are given by 



IT 



are the equilibrium distribution func- 



I(a, b, c, d) 



G(a, b, c, d) 



1 



+ » 



2 [(a 
1 



(6 + d) s 



sin 



2arctan 



2 [(a - c) 2 + (6 + d) 2 ' 

1 1 

2 [(a + c) 2 + (6 + d) 2 " 
1 1 



sin 



cos 



■ cos 



2arctan 
2arctan 
2arctan 



a + c 
a — c 



6 + d. 
a + c 



a — c 



6 + 



Notice that the function Giujr ,Tr ,0Jr , I~ 

y tR 1 IH' *L' *L 



(40) 

This resonant 



2[(a-c) 2 + (6 + rf) 2 

has a peak for cj^ 

behaviour is associated to the fact that the Higgs background is carrying a very 

low momentum (of order of the inverse of the bubble wall width L u ) and to the pos- 
sibility of absorption or emission of Higgs quanta by the propagating supersymmetric 
particles. The resonance can only take place when the left-handed stop and the right- 
handed stop do not differ too much in mass. By using the Uncertainty Principle, it is 
easy to understand that the width of this resonance is expected to be proportional to 
the thermalization rate of the particles giving rise to the baryon asymmetry. Within the 

3 To account for interactions with the surrounding plasma in the unbroken phase, particles must 
be substituted by quasiparticles, dressed propagators are to be adopted. This means that self-energy 
corrections ate one and twp-loops to the particle propagators modify the dispersion relations, which 
becomes u> 2 = k 2 + m 2 (T) with m(T) the plasma mass, and introduce a finite width T. 
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MSSM, however, it is assumed that the left-handed stop mass mj is much larger than 
the temperature T and the resonance can only happen for momenta larger than mqr . 
Such configurations are exponentially suppressed and do not give any relevant contri- 
bution to the CP-violating source. Nonertheless, if the electroweak phase transition is 
enhanced by the presence of some new degrees of freedom beyond the ones contained 
in the MSSM, e.g. some extra standard model gauge singlets, the resonance effects in 
the stop sector might be relevant. 

-The inelastic interactions of the right-handed stop number- 
Let us now investigate the other terms in Eqs. (30) and (31) leading to particle 
number changing interactions. The inelastic scattering processes which change the 
right-handed stop number are the ones induced by the top Yukawa coupling and by the 
axial top number violation inside the Higgs bubble wall (besides the the ones arising 
from strong sphaleron transitions). For sake of simplicity, let us focus on the top Yukawa 
interactions -other interactions may be treated similarly. If we assume that the left- 
handed stop is heavier than the temperature of the plasma and is therefore decoupled 
from the thermal bath, these interactions read 

C = htH$ P L tt* R + h.c, (41) 

where Pl denotes the left-handed chirality projector operator and H® and t are the 
Higgsino and top quark fields, respectively. 

The right-hand side of Eq. (|28[) contains now generic expressions like 

<S£S*Gfl + E£ 5G 7 - Gfi <5£f nt - 5G 7 E£ . (42) 

*R tR t R t R tjt t R tfj tR v 7 

Since the contribution to E~* induced by the top Yukawa quark coupling is symbolically 

t R 

given by hfG^ G tL , the expression fl42| ) may be written as 

sr/~i /~~t0 /~~t0 I s-iO r/~r /~~t0 I /~<0 /~iQ 

- G h SG H° G t~ G TR G % SGt L- SG t R G h 2 G t L - ( 43 ) 

We now show that -under a number of reasonable approximations- these terms all 
together lead to the familiar "scattering term" proportional to (nj R — n tL + ti^ j in the 
right-hand side of the diffusion equation for the right-handed stop number. 
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-Getting the diffusion equation- 



To progress further we need to introduce reasonable approximations based on physi- 
cal considerations. The main phenomena occuring in the plasma when the Higgs bubble 
wall is passing through a certain region are characterized by three different time (or 
length) scales: the quantum (microscopic) time tq, the statistical (macroscopic) time 
Ts and the time set by the bubble wall thickness and velocity r w ~ L^/v^. The first 
measures the range of radiative corrections to the Compton wavelength of particles and 
the second measures the range of interactions among particles. In the more familiar 
classic kinetic interaction theory, a similar distinction is made between the scattering 
length and the mean free path. It is a reasonable assumption that the statistical scale t$ 
is much larger than the quantum one. Furthermore, in the thick wall limit, the follow- 
ing hierarchy holds: ^> t$ ^> tq. This observation allows us to recast the quantum 
field theoretical problem into a much simpler kinetic theory and, under some further 
assumptions, to rederive the classical diffusion equations -using well-known techniques 



from relativistic many-body theory [16]. Suppose we separate the space-time "cells" 
whose characteristic length scale is intermediate between the statistical and the quan- 
tum scales. As the correlation between cells will be negligible by construction, the only 
interesting case is when the two arguments of a given propagator or self-energy lie in 
the same cell. In the interior of a single cell, relaxation phenomena are negligible. More 
concretely, the propagators may be Fourier transformed over a cell. Relaxation phe- 
nomena become apparent as we move from cell to cell. This picture breaks down if, 
for instance, the major contribution to CP-violating sources come from long wavelength 
particles for which the mean free path becomes comparable to the Compton wavelength. 
We will return to this point later. 

In mathematical language, if x and y are the arguments of a propagator or a self- 
energy, we Fourier transform with respect to the r = x—y, while the absolute coordinate 
X = |(x + y) serves as a cell label. The Fourier transform reads 

J wr^' ri " J,G (^) -I Wf^" G(x ' v) ' (44) 

where G(X,p) is called the Wigner transform of G(x,y). Similar formulae hold for the 
self-energies. 
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Let us now have a closest look at the third term in Eq. (^3]). We can rewrite it as 



T^ r (X,x 3 )SG Ir (x 3 ,X) 



iq-(X-x 3 ) 



x 5G 




e ik<S 3 -X)yO ( X + %3 rp . ; 



T, tg, fc 



(45) 



The integrand will be appreciably different from zero only when X belongs to the same 
cell as x 3 , that is (X + x 3 )/2 ~ X. Inserting the above expression into Eq. fl28|), we get 



dU 



d 3 k 



lS(X,q,T-t 3 )SGr(X,k,t 3 -T 



(46) 



( 27r )3 t« 

A similar expression is obtained for the last term in ( |43|) JG^E- . 

Let us now adopt the quasiparticle Ansatz for the shift in the right-handed stop 
Green function 



8G* (X, k) 



VK 



tR 



5Gf(X,k) = — 



l + 5f TR (X,k))6(k"-u^ R ) 

+ 5f TR (X,-k)5(k + u^ R ) , 
^(X,^ -^) 
+ (l + 54(X,-fc))5(A; + ^); , 



(47) 



where u~ = k 2 + to| (X, T) and the plasma mass squared rn~ (X, T) is a function of 

tR t R 1 t R 

the spacetime coordinate X as well as of the temperature T. This Ansatz is the logical 
generalization of the Green function for free fields and incorporate the renormaliza- 
tion effects of the non-trivial spectral density in the plasma. In this sense the Wigner 
functions 5fj R and 8fj are the quantum kinetic extension of the classical particls phase- 
space distributions for the right-handed stop and its antistate. Notice,though, that this 



Ansatz does not solve the QBE unless one neglects off-shell terms |23j. This procedure 
is however justified if t w ^> t s ^> tq. 

From this Ansatz we can construct the equation for the right-handed stop number 
asymmetry nj by inserting ( pTFD into the right-hand side of Eq. (^). After some algebra 
and exploiting the symmetry properties of the Green functions, we find that the third 
and the last term of Eq. (E3) lead to 



dn 



tR 



dT 



f V ■] 



tR 



oc 



X COS 



dti / du 



d 3 k 



(2tt) 3 tR 



*AJ(M S ), 



(48 
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where we have made use of the fact that the equilibrium self-energy E° does not depend 
on the "cell" label X and 

M)-Ej< L,k 



r~ k) = ^ V ! > ±±L± . (49) 



ILOr 

tR 

The kinetic Eq. (|48|) has an obvious interpretation in terms of gain and loss processes 
which change the number density asymmetry of the right-handed stop. This is because 
the rate of change of the particle number asymmetry in the plasma is related to the 
imaginary part of the retarded self-energy function, Eq. ([49]) . However the kinetic 
equation is non-Markovian and the retarded time integral and the cosine function replace 
the more familiar energy conserving function present in the classical kinetic equation. 
Under which conditions can we obtain the more familiar classical Boltzmann equation? 
Let us suppose first that the relaxation time scale for nj R is much longer than the time 
scale of the non-local kernel. Under this assumption, 5fj R (t 3 ) — 5fj R (t 3 ) can be replaced 
by 8fj (T) — 5 fa (T) and taken outside the time integral. This leads to a Markovian 
description. A further approximation is obtained taking the upper limit of the time 
integral to T — > oo. The cosine function becomes an energy conserving delta function 
and we obtain 

dnr -» r d 3 k 

(2tt) 

This result is the familiar relationship between the relexation time of the particle distri- 



h v 4 « - / 7^3 r S k' K - 6 k) & n ( 5 °) 



bution rs° and the damping rate, which is determined by the contribution to imaginary 
part of the self energy on shell from the top Yukawa coupling, Er = —i ( S~ 5> — S-^ ) lur . 

We can repeat the same procedure for all the other terms present in the expression 
( |43[) . Making the further approximation that the relaxation rate induced by the top 
Yukawa coupling is the same for the left-handed top, the right-handed stop and the 
higgsinoQ- r top - and that, in the high temperature limit, it does depend on the three- 
momentum only weakly , we can take it outside of the momentum integral. Moreover, 
if we follow the spirit of the usual derivation of Fick's law, we end up with the familiar 
diffusion equation for the right-handed stop number asymmetry 

r - D T R ^ 2n T R = "I** (nj R - n tL + n^) +■% + ••-, (51) 



dT 



4 This approximation may turn out to be very rough for some choice of the parameters because the 
imaginary part of the two-point self-energy of different particles depend sensitively on their dispersion 
relations in the plasma (23, E6| and, especially for Majorana fermions, the latter are highly nontrivial 
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where 

d 3 k 



717 



{ 5 h R ~ *4) » ( 52 ) 



'* J (2tt) 3 

Dj is the right-handed stop diffusion coefficient and the dots mean that one should 
also include the axial top number violation processes as well as the strong sphaleron 
interactions. 

To assess whether the Markovian and the Boltzmann approximations are reliable 
in the framework of supersymmetric electroweak baryogenesis one has to understand 
the different time scales. The time scale of the kernel is of the order of uj- 1 . Since 
Tj r <C uj r , the relaxation time scale of the right-handed stop population is longer than 
the range of the kernel and the Markovian approximation is warrented. Furthermore, 
since the time scale set by the bubble wall profile ~ L w /v UJ and the diffusion time scale 
~ Dj r /v^ are larger than the relaxation time, the Boltzmann/diffusion equation is also 
warrentedQ 

6. The QBE for the Higgsino number 

The Higgs fermion current associated with neutral and charged Higgsinos can be 
written as 



Jt = Hj^H (53) 

H 



where H is the Dirac spinor 



*=if;) (54) 

and H 2 = H 2 (Ht)i -^1 = -^1 (Hi) f° r neutral (charged) Higgsinos. The processes 
in the plasma which change the Higgsino number are the ones induced by the top 
Yukawa coupling and by interactions with the Higgs profile. The interactions among 
the charginos and the charged Higgsinos which are responsible for the CP-violating 
source in the diffusion equation for the Higgs fermion number read 

C = -g 2 [vi(x)P L + e ie »v 2 {x)P B \ w} + h.c, (55) 

where 9^ is the phase of the /^-parameter. Analogously, the interactions among the 
Bino, the H^-ino and the neutral Higgsinos are 

£ = -I jfo [ Vl ( x )p L + e ^v 2 (x)P R ] (g 2 W 3 - 9l B)} + h.c. (56) 



5 On the other hand, for soft scales and for values of the right-handed SUSY breaking mass term such 
that the plasma mass m—^ is small compared to the temperature, the relaxation time scale becomes 
comparable to the time scale of the kernel and the Markovian approximation breaks down. 
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We can compute the QBE for the Higgsino number repeating step by step the procedure 
adopted in the previous section. To compute the source for the Higgs fermion number 
we again perform a "Higgs insertion expansion" around the symmetric phase. At the 
lowest level of perturbation, the interactions of the charged Higgsino induce a contribu- 
tion to the self-energy of the form (and analogously for the other component <5S CP 



H 



S^ix, y) = g% P {x, y)P L G%<{x, y)P L + g*p(x, y)P R G^{x, y)P R 



W 



w 



(57) 



where 



9cp( x ,y) 



g 2 2 v 1 (x)v 2 (y)e i6 » 
giv 1 (y)v 2 (x)e iB ' 1 . 



(5f 



Again, we have approximated the exact Green function of winos by the equilib- 
rium Green function in the unbroken phase G°~ . Similar formulae hold for the neutral 

w 

Higgsinos. 

Analogously to the case of right-handed stops, the dispersion relations of charginos 
and neutralinos are changed by high temperature corrections [^J . Even though fermionic 
dispersion relations are highly nontrivial, especially when dealing with Majorana fermions 
relatively simple expressions for the equilibrium fermionic spectral functions may 



be given in the limit in which the damping rate is smaller than the typical self-energy 
of the fermionic excitation |17| . The computation of the source goes along the same 
lines of the calculation done in the previous section and it is easy to show that the 



CP- violating source [15 



H 



- / d 3 x 3 



dt* Tr 



+ G°f(X,x 3 )5^ h 



5Xf><(X,x 3 )GY(x 3 ,X) 



""^(ss,*) 

containes in the integrand the following function 



i0,>, 



(V ^X3,X) 

(59) 



g CP (X, x 3 )+g CP (X, x 3 )-g CP (x 3} X)-g CP (x 3 , X) = 2ism9 / ,[v 2 (X)v l (x 3 ) - fi(X)u 2 (x 3 )] 

(60) 

which vanishes if Im(/i) = and if the tan/?(x) is a constant along the Higgs profile. 
Performing the "Higgs derivative expansion", we finally get 



3M 2 g\ Tf + M x g{ X| 



2 tB 



(61) 
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where 



dk- 



k 2 



2lT 2 U7iUJ 



2Re(f\ 





w 



r^) + (l-2Re(/|))/(^,r 



+ 



2(lm(f s )+Im(f~))G(u s ,r~., 



alld 4(H>) 



k 2 + |/i| 2 (Mf ) while /f ( ~ } 



(62) 
The 



1/ [exp (u s ^/T + i?H{w)l T ) + 1 
exact computation of the damping rate of charged and neutral Higgsinos will be 
given elsewhere [EB|, but it is expected to be of the order of 5 x 10~ 2 T. The Bino 
contribution may be obtained from the above expressions by replacing M2 by Mi. As 
for 5^- , the CP- violating source for the Higgs fermion number is enhanced if M2, Mi ~ /x 
and low momentum particles are transmitted over the distance L u . This means that the 
classical approximation is not entirely adequate to describe the quantum interference 
nature of CP-violation and only a quantum approach is suitable for the computation 
of the building up of the CP- violating sources |15| . 

To find the "collision terms" which change the higgsino number inside the Higgs 
bubble wall one has to repeat the analysis performed previously for the right-handed 
stop number. In particular, one has to adopt the quasiparticle Ansatz for the Higgsinos. 
The one for charged Higgsinos, for instance, reads 



SG> s (X,k) 



-^( / k+\ fl \)[fl + 6f s (X,k))6(k° 
+ 5/~(X,-£)5(£; + ^; 



CO 



H 1 



in 



5G^(X,k) = U+M) 5f S (X,k)5(k»-u s ) 



H 



+ (l + 5f s {X,-k)) 6{k°+u s ) 



(63) 



where we have assumed that the mass term is well approximated by the bare mass term 
Inserting this Ansatz into the QBE for the fermionic degrees of freedom ( p5|) one 
can recover the usual Boltzmann/diffusion equation for the Higgsino number inside the 
Higgs bubble wall once the same approximations discussed in the subsection 5.2 have 
been made. 
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7. The baryon number asymmetry 

Once we have computed the CP- violating sources for the right-handed stop number and 
the Higgsino number and we have shown that in the thick wall limit the particle number 
changing interactions in the QBE's reduce to the more familiar Boltzmann/diffusion 
equations, we are ready to estimate the final baryon asymmetry produced during the 
electroweak phase transition. 

From now on we will closely follow the approach taken in ref. [|7|, || where the reader 
is referred to for more details. If the system is near thermal equilibrium and particles 
interact weakly, the particle number densities rij may be expressed as rii = kifiiT 2 /Q 
where /ij is the local chemical potential, and ki are statistical factors of order of 2 (1) for 
light bosons (fermions) in thermal equilibrium, and Boltzmann suppressed for particles 
heavier than T. 

The particle densities we need to include are the left-handed top doublet qi = 
(tL + the right-handed top quark t R , the Higgs particle H = (H®, H®, , H2), 
and the superpartners t R and H -remember that the left-handed stop and sbottom are 
supposed to be heavier than the temperature T. As usual, we shall assume that the 
supergauge interactions are in equilibrium. Under these assumptions the system may be 
described by the densities Q = qL, T = tn + tn and 7i = H + H . Ignoring the curvature 
of the bubble wall, any quantity becomes a function of the coordinate z = z 3 + v^z, 
the coordinate normal to the wall surface, where we assume the bubble wall is moving 
along the Z3-axis. 

Assuming that the rates of the interactions induced by the top Yukawa coupling r top 
and by the strong sphalerons r ss are so fast that Q/k q — H/k^ — T/kr = 0(l/T top ) 
and 2Q/k q — T/kr + 9(Q + T)/kb = 0(1/T SS ), one can find the equation governing the 
Higgs density 

vjH - DH" + TH - 7 = 0, (64) 

where the derivatives are now with respect to z, D is the effective diffusion constant, 
7 is an effective source term in the frame of the bubble wall and T is the effective 
decay constant 0. An analytical solution to Eq. (|^) satisfying the boundary condi- 
tions H(±oo) = may be found in the symmetric phase (defined by z < 0) using a 
z-independent effective diffusion constant and a step function for the effective decay 
rate T = T6(z). A more realistic form of T would interpolate smoothly between the 
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symmetric and the broken phase values. One can check, however, that the result is 
insensitive to the specific position of the step function inside the bubble wall ||. The 
values of D and r in (53) of course depend on the particular values of supersymmetric 



parameters. For the considered range we typically find D ~ 0.8 GeV -1 , V ~ 1.7 GeV 
§■ 

The tunneling processes from the symmetric phase to the true minimum in the first 
order phase transition of the Higgs field in the MSSM has been recently analyzed in 
including the leading two-loop effects. It was shown that the Higgs profile along the 
bubbles at the time when the latter are formed have a typical thickness ~ (20 — 30) /T 
f\. The total amount of the baryon asymmetry is proportional to A/3 -the change in 
the ratio of the Higgs vacuum expectation values (3 = v%/v\ from z = to inside the 
bubble wall. This quantity tends to zero for large values of rriA, and takes small values, 
of order 10 -2 for values of the pseudoscalar mass rriA = 150-200 GeV [2? |. 



The dependence of the final baryon asymmetry on 7 is much more delicate. The 
solution of Eq. (|64[) for z < is 

H{z) =Ae zv - /15 , (65) 



and for z > is 




1 



D(X + - A_) A) 

_J f 

D(A_ - A+) Jo 



du^(u)e- x+u j e A+z 
' du^(u)e' x - u )e x - z . (66) 



where 

2D y 1 

and 7(z) is the total CP- violating current resulting from the sum of the right-handed 

stop and Higgsino contributions. Imposing the continuity of 7i and Hi at the boundaries, 

6 In general, however, the value of L u when the bubbles are moving through the plasma with some 
velocity is different from the value at bubble nucleation. Indeed, the motion of the bubble wall is 
determined by two main factors, namely the pressure difference between inside and outside the bubble 
-leading to the expansion- and the friction force, proportional to v u , accounting for the collisions of 
the plasma particles off the wall. The equilibrium between these two forces imples a steady state with a 
final velocity v u . If bubbles are rather thick, thermodinamical conditions are established inside the wall 
and for the latter is no longer possible to loose energy by thermal dissipation. Under these conditions 
the bubble wall is accelerated until slightly out-of-equilibrium conditions and the friction forces are 
reestablished. 



24 



we find 



A = B + (l-±) =B-(±-l\ = =^- /_" "d«7(«)e-^. (68) 



From the form of the above equations one can see that CP-violating densities are non 
zero for a time t ~ D/v^ and the assumptions leading to the analytical form of W(z) 
-from Eq. (64)-are valid provided that the interaction rates r top and T ss are larger than 

v 2 jd e a- 

The equation governing the baryon asymmetry tib is given by |7| 

^Vb - v w n' B - 9(-z)n f T ws n L = 0, (69) 



where T ws = Qna%T is the weak sphaleron rate (k ~ 1) |8| (the correct value of « is 
at present subject of debate), and is the total number density of left-handed weak 
doublet fermions, rif = 3 is the number of families and we have assumed that the baryon 
asymmetry gets produced only in the symmetric phase. Expressing n^(z) in terms of 
the Higgs number density 

_ 9k q k r - 8k b k T - 5k b k q 
HL ~ k n (k b + 9k q + 9k T ) [ } 

and making use of Eqs. fl65l)-(|69D, we find that 



— = -9 {kij — 5 , 71 

s v^s 

where s = 2ir 2 g* s T 3 /45 is the entropy density (g* s being the effective number of rela- 
tivistic degrees of freedom) and g(ki) is a numerical coefficient depending upon the light 
degrees of freedom present in the thermal bath. 

We see that the final baryon asymmetry depends sensitively on the parameter A, 
that is on the integral of the source. In previous analysis [0, the CP-violating sources 
have been inserted into the diffusion equations by hand only after the CP-violating 
currents have been defined and computed. More specifically, CP-violating sources S 
associated to a generic charge density j° have been constructed from the current j M by 
the definition S = doj° J7], |]. This procedure has introduced an unacceptable degree 
of arbitrariness -different definitions of CP-violating currents j M have been adopted in 
the literature 0, § and have lead to different numerical results for the final baryon 
asymmetry. This is because the corresponding sources differ in terms of number of 
derivatives of the Higgs bubble wall profile and, therefore, in terms of different powers 
of the bubble wall velocity v u and bubble wall width L w . 
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The basic merit of the CTP formalism is to provide us with a rigorous and self- 
consistent definition of the CP-violating sources within the quantum Boltzmann equa- 



tions [pTo| . The computation of the CP- violating currents for the right-handed stop and 
higgsino local densities was performed in || by means of the CTP formalism. Since 
currents were proportional to f3 in the thick bubble wall limit, i.e. proportional to the 
first time derivative of the the Higgs profile, sources turned out to be proportional to the 
second time derivative of the Higgs profile |8] . However, the derivation of the transport 
QBE's allows us to define the CP-violating sources uniquely, since they are automatically 
incorporated in the right-handed side of the QBE's fl5|]. This self-consistent procedure 
indicate that the sources in the quantum diffusion equations are proportional to the 
first time derivative of the Higgs configuration. A comparison between the sources Sj r , 
see Eq. (|38D , and Sg, see Eq. (|S"T|) , obtained in the present work and the currents j° 
given in ref. || indicate that they may be related as 

j°(T) 

S{T) ~ ±-LJ. (72) 

and it may be interpreted as the time derivative of the current density accumulated at 
the time T after the wall has deposited at a given specific point the current density j° 
each interval r 

S(T) ~d J T dt 3 -^. (73) 

Here r = T _1 is the thermalization time of the right-handed stops and higgsinos, re- 
spectively. The integral over time is peculiar of the quantum approach and it induces 
memory effects. 

The parameter A computed from the sources fl38|) and (|61|) and rigorously derived 
from the QBE's turns out to be(] 

2f(ki)T~ r°° n / \ \ 2/(A;i)r~ 
D \ + Jo H D X + 

I = rduv\u) d -^e-^^ rduv\u) d ^ (74) 
Jo du Jo du 

where j ~ is the Higgsino current computed in || , is the higgsino damping rate and 
f{ki) is a coefficient depending upon the number of degrees of freedom present in the 
thermal bath. 



7 For large values of the left-handed stop mass the stop contribution to the baryon asymmetry is 
strongly suppressed compared to the chargino and neutralino ones and the numerical values of the 
baryon asymmetry depend linearly on the phase of the Higgsino mass parameter. 
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A comparison between this parameter A and the one obtained in M indicates that 
memory effects introduce an enhancement factor in the final baryon asymmetry 

" (75) 



X+v^ \5 x 10- 2 7y \v L _ 

with respect to the results given in || 0. Our results indicate that memory effects play a 
crucial role in relaxing the stringent lower bound on the values of the CP- violating phases 
obtained in previous analysis f|, |10|| . In particular, the observed baryon asymmetry 
?ib/s~4x 10~ n may be originated by Higgsinos if the CP- violating phase of the 
/i-parameter is 



v. 



|sin(0,)|>lO-^^j. (76) 
This lower bound is two orders of magnitude weaker than the one obtained in previous 



works H |T0[- These small values of the phase M are consistent with the constraints 



from the electric dipole moment of the neutron even if the squarks of the first and second 
generation have masses of the order of 100 GeV. In the limit of thick bubble walls, the 



CP- violating sources are also characterized by resonance effects || pl| when the particles 
involved in the construction of the source are degenerate in mass. The resonance is 
manifest in the function G defined in fl40|). The interpretation of the resonance is 
rather straightforward if we think in terms of scatterings of the quasiparticles off the 
advancing low momentum bubble wall configuration. For momenta of order of the 
critical temperature, the scattering is more efficient when the Higgsinos and gauginos 
are nearly degenerate in mass, \x ~ M 2 . A similar effect has been found in ref. [ ]TTJ 
where the system was studied in the classical limit. These classical treatments somehow 
obscure the origin of the CP-violating effects as resulting from quantum interference and 
the origin of the resonance. Moreover, memory effects are not present in the classical 
approach which, therefore, understimates the value of the final baryon asymmetry. 



8 The integral I has been computed including two-loop effects in ref. J27J and results to be / ~ 10~ 2 
for uiA — 150-200 GeV. The same value (in absolute value) of the integral was obtained in ref. § 
by making use of a reasonable Ansatz for the Higgs profile. The reason is that this result was based 
on the temperature Tq defined to be the temperature at which the Higgs potential becomes flat at the 
origin. This temperature is lower than the temperature T c at which the tunneling process occurs and 
the numerical result of ref. 27 indicates that the integral I is a decreasing function of the temperature. 
Since the results of ref. || were based on the one-loop effective potential, the result in |^7| is accidentally 
very close to the numerical result of ||. Therefore, the numerical difference between the value of the 
final baryon asymmetry obtained here and the one evaluated in ref. |^| reduces to the expression j75|). 
We thank M. Quiros for discussions on this point. 
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8. Conclusions 

In conclusion, we have reanalyzed the computation of the baryon asymmetry gen- 
erated at the electroweak phase transition in the Minimal Supersymmetric Standard 
Model. We have argued that the CTP formalism is the natural guide towards a rig- 
orous and self-consistent definition of the CP-violating sources within the quantum 
Boltzmann equations. The QBE's manifest -in the scattering terms -non-Markovian 
features. In the limit of thick bubble wall, however, a number of approximations may 
be performed so that the QBE's get reduced to the more familiar diffusion/Boltzmann 
equations. Strong memory effects enhance the strength of the CP-violating sources by 
at least two orders of magnitude. The baryon asymmetry is mainly generated by Hig- 
gsinos, provided that Higgsinos and gauginos are not much heavier than the electroweak 
critical temperature and the phase M is larger than 1CT 3 . This value is further reduced 
for v u < 0.1. It is intriguing that these small values of the phases are perfectly consis- 
tent with the constraints from the electric dipole moment of the neutron and squarks 
of the first and second generation as light as ~ 100 GeV may be tolerated. 
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